Classical Phase Fluctuations in High Temperature Superconductors 
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Phase fluctuations of the superconducting order parameter play a larger role in the cuprates than 
in conventional BCS superconductors because of the low superfluid density ps of a doped insulator. 
In this paper, we analyze an XY model of classical phase fluctuations in the high temperature 
superconductors using a low-temperature expansion and Monte Carlo simulations. In agreement 
with experiment, the value of ps at temperature T = is a quite robust predictor of Tc, and 
the evolution of ps with T, including its T-linear behavior at low temperature, is insensitive to 
microscopic details. 

Two classes of thermal excitations are responsible for 
disordering the ground state of a superconductor: fluctu- 
ations of the amplitude and phase of the complex order 
parameter. A consensus has not yet been reached on the 
relative importance of the two in the high temperature 
superconductors, since both are anomalous. The low su- 
perfluid density (phase stiffness) of the doped insulator 
impHes that phase fluctuations play an unusually large 
role. [|l[ Yet the nodes in the d-wave gap function sup- 
port more quasiparticle (amplitude) excitations at low 
temperatures than in a clean s-wave superconductor. 

This paper is concerned with an analytical and numer- 
ical study of the thermal evolution of the in-plane helicity 
modulus, 7||(r), of an anisotropic quasi two-dimensional 
classical XY model of phase fluctuations in a high tem- 
perature superconductor, pi We neglect quasiparticle 
fluctuations because the nodal quasiparticles are excita- 
tions of the insulating state; very little charge transport 
is associated with them, and their contribution to the 
superfluid density should be proportional to a positive 
power of the (small) doping concentration. B We also 
neglect collective amplitude fluctuations associated with 
the quantum dynamics of the phase since, with sufficient 
screening, the phase fluctuations are predominantly clas- 
sical down to quite low temperature. H 

The calculations focus on the scaled curve, 7|| {T)/^\\ (0) 
vs. T/Tc, and the value of the dimensionless ratios 
Ai = r,/7||(0) and A2 = T,7[|(0)/7||(0). (Here 7(1 (T) = 
d'y\\{T) / dT .) These nonuniversal quantities turn out to 
be rather insensitive to microscopic details of the model, 
such as the strength of the interplane coupling and the 
exact short-distance nature of the interactions, as shown 
in Fig. 1 and in the tables. Figure 1 also shows that the 
model results agree well with experiment, H when the 
helicity modulus of the model is related to the in-plane 
superfluid density, ps as determined by 



FIG. 1. Superfluid density vs. temperature, scaled by the 
zero-temperature superfluid density and by Tc, respectively. 
Experimental data on YBCO is depicted by the black line, 
and is taken from Hardy et al. 0] (The data are essentially the 
same for a range of doping concentration.) Our Monte Carlo 
results for system size 10 x 10 x 10 are the unfilled symbols. 
Calculations are for two planes per unit cell, with coupling 
J|l = 1 within each plane, and J^ and J^ between alternate 
planes, as defined in Eq. (3). Monte Carlo points above Tc 
are nonzero due to finite size effects. 
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where aj_ is the spacing between planes and \ab is the 
London penetration depth within the Cu02 planes. 



There is strong empirical evidence that classical phase 
fluctuations determine much of the important physics in 
the superconducting state of the high Tc superconduc- 
tors, and also some properties of the normal state, espe- 
cially in underdoped materials. 0| Most notably, T^ in- 
creases roughly linearly with the zero temperature super- 
fluid density, [Ai = rc/7||(0) ^ 1) whereas the char- 
acteristic energy scale for pairing, 1^0/"^, is both quan- 
titatively large compared to ksTc and decreases as the 
doping increases. Furthermore, ARPES and other mea- 
surements of the superconducting gap reveal that pair 
formation occurs at a crossover temperature well above 
Tc- [|[-|0| It is important to note that A2 = T^i^^ {0)h\\ (0) 
is roughly constant for various materials and doping con- 
centrations. This implies that the fluctuations predom- 
inantly responsible for the T-linear dependence of the 
superfluid density at low T are also responsible for the 
ultimate destruction of the superconducting state at Tc. 
As first pointed out by Roddick and Stroud fl^], this 



behavior is characteristic of classical phase fluctuations. 

The following arguments have been made against this 
interpretation of the data: 1) The non-universal ratio 
Ai = rc/7||(0) should ||l^ theoretically lie in the range 
4-8, rather than in the experimentally observed range of 
0.5 - 1. In particular, it has been argued that Ai (x n 
in multilayer materials, where n is the number of layers 
per unit cell. 2) For weakly coupled layers, a phase only 
model would [Of yield a Ps{T) curve that looks like a 
rounded Berezinskii-Kosterlitz-Thouless [BKT) discon- 
tinuity, unlike what is seen in experiments. 3) Quan- 
tum effects suppress classical phase fluctuations Q for 
temperatures below the plasma frequency. 4) If pairing 
occurs in a substantial range of temperatures above Tc, 
the effects of fluctuation superconductivity should be ob- 
served, contrary to experiment. [[l2| 

As can be seen from the figure and the tables, our 
present results conclusively show that the first two as- 
sertions are incorrect; the classical XY model is quan- 
titatively consistent with experiments. The third sug- 
gestion has been previously shown to be incorrect [g, 
due to screening by the substantial background normal 
fluid. Specifically, in a two-fluid model of a superconduc- 
tor, the classical model is reliable down to a classical to 
quantum crossover temperature which can be well below 
Tc'- Tciass oc Tc/<jn, where a^ is an average of the op- 
tical conductivity of the normal component in units of 
the quantum of conductance. The fourth point overlooks 
the fact that fluctuation superconductivity is only signif- 
icant close to Tc where the correlation length is long. In 
conventional superconductors the observed fluctuations 
involve amplitude and phase and are Gaussian, while for 
the high Tc superconductors, true critical fluctuations in 
the XY universality class are detected in a remarkably 
broad range of temperatures. ||l5|| 

A classical XY model on a tetragonal lattice will be 
used to study the effects of phase fluctuations in a quasi 
two-dimensional superconductor at wavelengths that are 
long enough for amplitude fiuctuations to be unimpor- 
tant. The in-plane unit cell area an does not enter into 
the evaluation of Tc or the temperature dependence of 
7||(r). Here a|| is a short-distance cutoff which will be 
discussed at the end of the paper. In general, the in- 
teraction energy, V, depends on the phase difference, 
6ij = 9i ~ 9j, between nearest-neighbor sites < i,i >. 
Because of gauge invariance and time reversal symmetry, 
V can be expanded in a cosine series. 



^(^u) = X! ^" cosindi. 
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The first harmonic, cos(6'), corresponds to the transfer 
of one pair of electrons between neighboring cells; each 
successive harmonic transfers a higher number of pairs. 
We keep only the first two terms in the cosine series for 
couplings within a plane, and the first cosine term for the 
weaker Josephson coupling between planes. 



H = -J|| Y^ {cos(%) +Scos{2eij)} 

- ^ {jfcos(0fcO}, (3) 

where the first sum is over nearest neighbor sites within 
each plane, and the second sum is over nearest neigh- 
boring planes. The coupling, J|| , will be assumed to be 
isotropic within each plane and the same for every plane, 
but the coupling between planes, J'f , is different for crys- 
tallographically distinct pairs of neighboring planes. It 
will be assumed that Jii, J±, and S are positive, since 
there is no reason to expect any frustration in the prob- 
lem, [|l6| and that S < 0.25, since for S > 0.25 there 
is a secondary minimum in the potential for 9ij = tt, 
which is probably unphysical. The sensitivity of various 
computed quantities to variations in S in this range is a 
measure of the importance of "microscopic details." 

It follows from simple and general considerations that 
most features of the thermal evolution of the superfluid 
density of YBa2Cu307_5 shown in Fig. 1 are reproduced 
by such a model. The critical phenomena are in the 
same universality class as the classical 3D XY model, 
which is consistent with the observed behavior |15[| in 
YBa2Cu307_5 within about 10% of Tc. Furthermore, 
the helicity modulus is linear in the temperature, as ob- 
served in YBa2Cu307_5. Indeed, using linear spin-wave 
theory, it is straightforward to obtain the first terms in 
the low-temperature series for 7: 



711 (T) = J|| (1 + 4(5) 
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where (a; is a numerical integral which varies from a = \ 
in the 2d limit (J_i -^ 0), to o; = 2/3 for J^ -^ 7||(0) = 
J|| (1 + 4(5) for one plane per unit cell. 

A more quantitative comparison between the classi- 
cal model and experimental data can be undertaken by 
studying various dimensionless ratios, particularly Ai = 
Tc/7||(0) and A2 = Tc7[|(0)/7||(0). Here, Tc is computed 
numerically by means of the Binder parameter |lj| for 
systems of size up to 24 x 24 x 24. Errors in Tc are 
limited by the resolution with which the Binder crossing 
point is computed in each case. The quantities 7|| (0) and 
7|l (0) are obtained from Eq. (^) . 

Experimentally A2 ~ |, and Ai is in the range 0.6-1.3 
for underdoped and optimally doped materials. Tables 1 
and 2 show the ratios Ai and A2 for various choices of 
parameters in the classical XY model. Note that A2 ^ ^ 
for (5 not too small, whereas A2 is about a factor of two 
smaller for S — 0. The shape of the Ps{T) vs. T curves, 
as quantified by v42, is remarkably robust, especially if 
we compare the cases of (5 = 0.1 to 5 — 0.25. The ratio 
Ai is a little more sensitive to the value of (5, but it is 
comfortably in the experimental range for 5 between 0.1 
and 0.25, and only slightly larger for (5 = 0. The relative 



insensitivity of Ai to 5 and to the details of the interplane 
couplings, J_L, demonstrates that when classical phase 
fluctuations govern the physics, Ps(0) is a quantitatively 
good predictor of Tc- 
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0.1 


0.25 





0.1 


0.25 





0.1 


0.25 


Ai 


0.89 


0.72 


0.6 


1.1 


.828 


.625 


1.324 


.986 


.73 


A2 


.22 


.335 


.38 


.27 


.381 


.388 


.3066 


.432 


.437 


Tc 


.89 


1.01 


1.2 


1.1 


1.16 


1.25 


1.324 


1.38 


1.46 



Table 1: Single Layer The dimensionless ratios Ai — 
Tc/-f\\{0) and A2 = T'c7||(0)/7||(0) which characterize the 
superfluid density vs. temperature for systems with one 
layer per unit cell. Jj_ and T^ are quoted in units of Jy . 
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0.1 


0.25 





0.1 


0.25 


Ai 


1.38 


1.03 


0.78 


1.13 


0.836 


0.645 


A2 


.279 


.402 


.426 


.271 


.376 


.394 


J± 


1 : .01 


1 : .01 


1 : .01 


.1 : .01 


.1 : .01 


.1 : .01 
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0.25 





0.1 


0.25 


Ai 


1.52 


1.12 


0.83 


1.2 


0.907 


0.675 


A2 


.306 


.437 


.452 


.29 


.407 


.411 


J± 


1 : .1 


1 : .1 


1 : .1 


.1 : .1 


.1 : .1 


.1 : .1 


d 





0.1 


0.25 





0.1 


0.25 


Ai 


1.6975 


1.252 


0.916 


1.324 


0.986 


0.73 


A2 


.3318 


.4772 


.4906 


.3066 


.432 


.437 



Table 2: Double Layer The dimensionless ratios Ai = 
Tc/j\\ (0) and A2 = Tc'y'n (0)/7|| (0), which characterize the 
superfluid density vs. temperature curve, for systems 
with two layers per unit cell; the two values of J± are 
between planes within a unit cell and between unit cells. 

Hardy et al. have demonstrated ||6| that when 
Ps{T)/ps{0) is plotted vs. T/Tc, for various dopant con- 
centrations in YBa2Cu307_5, the data collapse approx- 
imately onto one curve. Since Tc oc Ps(0), this amounts 
to scaling out Tc for both axes. Thus the unique shape 
of the normalized ps vs. T curve may be attributed to 
the existence of a single energy scale, the transition tem- 
perature. This is inconsistent with scenarios that require 
the superfluid density to be depleted by quasiparticle ex- 
citations as the temperature is raised, since the energy 
scale is set by Aq, and Aq/Tc is not roughly constant, 
but varies by a factor of 2 in the range of doping studied 
by Hardy et al. 

We have also used Monte Carlo calculations to evalu- 
ate the superfluid density for < T < Tc. The results are 
scaled as in Hardy et al. M and compared with their data 
in Fig. 1. As anticipated, for S not too near zero, the 
model is insensitive to the value of S, and in remarkable 
agreement with experiment, considering that no param- 
eters have been tuned. Quasiparticle excitations near 



the nodes of the d-wave gap may give a T-linear contri- 
bution to the superfluid density, but our results suggest 
that this contribution is quantitatively small. This, in 
turn, supports the idea that there is little charge trans- 
port associated with the quasiparticles. 

A much-discussed feature of the systematics of Tc in 
the high temperature superconductors pSJ is the ob- 
served increase of Tc within each family of materials with 
the number of planes per unit cell, n. Within the clas- 
sical phase model, the fact that phase fluctuations lead 
to a particularly large suppression of Tc below its mean- 
field value (see Table 3) leads to an increased sensitivity 
to even weak couplings in the third direction. This pro- 
duces a strong increase of Tc with n, although possibly 
not quite as strong as observed experimentally. However, 
it should be noted that other things may change with n; 
for example, in a three-plane material, the central plane 
may have a different hole concentration than the others. 
At the same time, since the dimensionless constant a in 
Eq. (^ is only weakly dependent on J±, A2 is an increas- 
ing function of n. This is to be expected as small values 
of A2 are characteristic of low dimensional phase fluctu- 
ations. Table 3 shows the variation of Tc with number of 
planes per unit cell for weak coupling between planes. 
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1 


2 


3 


4 


00 


dp(0) 
dT 


.2472 


.2384 


.2365 


.2348 


.2315 


Tc 


1.09 


1.20 


1.24 


1.26 


1.324 


Tmf 


1.111 


1.287 


1.334 


1.361 


1.394 
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0.1 


0.1 


0.1 


0.1 


0.1 


n 


1 


2 


3 


4 


00 


dp{0) 
dT 


.4604 


.4486 


.4450 


.4427 


.4379 


Tc 


1.16 


1.27 


1.30 


1.32 


1.38 


S 


0.25 


0.25 


0.25 


0.25 


0.25 


n 


1 


2 


3 


4 


CX) 


dpW 

dT 


.6211 


.6092 


.6055 


.6032 


.5982 


Tc 


1.25 


1.35 


1.39 


1.41 


1.46 



Table 3: Variations as a function of the number of 
planes per unit cell: The coupling between planes within 
the unit cell is Jj^ = 0.1, and between planes of different 
unit cells is J'j^ = 0.01. Tmf is the interplane mean-field 
estimate of Tc obtained as described in the text. 

Mean field theory is a standard method of estimating 
the effects of weak higher-dimensional couplings on the 
critical temperature of quasi one or two dimensional sys- 
tems. For instance, for one plane per unit cell (n — 1) 
this approach leads to an implicit equation for the three 
dimensional Tr 



X2diTMF)'2Jl_ = 1, 



(5) 



where X2d(T) is the susceptibility of an isolated plane. 
For the case 5 = 0, we have computed the interplane 



mean-field transition temperature, which is also pre- 
sented in Table 3, using the Monte Carlo results of Gupta 
and Baillie [^ for X2d{T). This mean- field theory be- 
comes exact in the limit J^ ^ 0, and always gives an 
upper bound to Tc- 

Phase fluctuations should also have detectable effects 
on other equilibrium properties, such as the specific heat, 
the diamagnetic susceptibility, and j± . In contrast to 7|| , 
these quantities depend on ay. Classical systems have 
nonzero specific heat at T = 0, in violation of the third 
law of thermodynamics. The classical XY model at tem- 
peratures T -^ Tc has a specific heat per unit area in a 
Cu02 plane equal to C = ks/^ai. The specific heat ||2^] 
at T = 2K of good crystals of optimally doped YBCO is 
roughly 5 x 10"'* /cs per planar copper; if we assumed that 
all of this specific heat were due to phase fluctuations, it 
would imply aii — 32 lattice constants. 

In the classical XY model a|| ~ r^, where r.y is the 
radius of a vortex core. Recent //SR measurements ||20|] 
have found that r„ grows substantially at low fields, and 
tends to a zero field value which is on the order of lOOA 
(26 lattice constants) and which is only weakly tempera- 
ture dependent nearly up to Tc- This behavior is consis- 
tent with other measurements, especially STM in a field 
of 6T, |2^] which have produced significantly smaller es- 
timates of the core radius. Thus, if we estimate ay using 
the /kSR data, it is consistent to attribute a large fraction 
of the low-temperature specific heat to classical phase 
fluctuations. The contribution of critical fluctuations to 
the specific heat near Tc may also be dominated by clas- 
sical phase fluctuations, but a quantitative comparison 
of the theoretically expected (non- universal) critical am- 
plitudes with experiment is not straightforward. 

Experimentally, p±/m* is much more weakly temper- 
ature dependent than p\\/m* at low temperatures. We 
intend in the near future to analyze the implications of 
this observation for the phase fluctuation model. Specif- 
ically, since 7_l/7|| = pj_/p||(a||/a^)^, the anisotropy of 
the superfluid density may also contain valuable informa- 
tion concerning the physics of aii . 

Finally, we address the remarkable measurements of 
the frequency dependent superfluid density in BSCCO of 
Corson et al. [E3[ Without making any explicit assump- 
tions concerning the dynamics, we can interpret these re- 
sults in terms of a finite size scaling hypothesis, in which 
we associate a length scale, L{uj), with the finite mea- 
surement frequency, and 



7(r,L)^L-T7(L/e(T)) 



(6) 



where ^(T) is the correlation length of the infinite system 
at temperature T. Since BSCCO is highly anisotropic, 
we follow Corson et al |23] in assuming that the finite fre- 
quency response is essentially two dimensional, in which 
case X = 0, and ^(r) is infinite for all T < Tbkt, the 
BKT transition temperature. This implies that 7 is ap- 
proximately frequency independent for T < Tbkt and 



vanishes exponentially as a function of L/^ at tempera- 
tures enough above Tkt that L > ^(T). Indeed, ^{T,L) 
was computed numerically for the two dimensional XY 
model by Schultka and Manousakis H and we have re- 
peated these calculations for anisotropic three dimen- 
sional models; the results confirm that our model nicely 



accounts for the observations of Corson et al. |23|, with 
the proviso that the measured finite frequency superfluid 
density is interpreted as a renormalized, rather than a 
"bare" , response function. 
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